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Abstract: We study Maxwell's equations in time domain in an anisotropic medium. 
The goal of the paper is to solve an inverse boundary value problem for anisotropies 
characterized by scalar impedance a. This means that the material is conformal, 
i.e., the electric permittivity e and magnetic permeability [l are tensors satisfying 
fj, = a 2 e. This condition is equivalent to a single propagation speed of waves with 
different polarizations which uniquely defines an underlying Riemannian structure. 
The analysis is based on an invariant formulation of the system of electrodynamics 
as a Dirac type first order system on a Riemannian 3— manifold with an additional 
structure of the wave impedance, (M,g,a), where g is the travel-time metric. We 
study the properties of this system in the first part of the paper. In the second 
part we consider the inverse problem, that is, the determination of (M,g,a) from 
measurements done only on an open part of the boundary and on a finite time 
interval. As an application, in the isotropic case with M C M 3 , we prove that the 
boundary data given only on an open part of the boundary determine uniquely the 
domain M and the coefficients e and \i. 

Keywords: Maxwell's equations, manifolds, inverse problems. 
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Introduction 



In this paper we study direct and inverse boundary value problems for Maxwell's 
equations on Euclidean domains in M 3 and on compact manifolds. In a bounded 
smooth domain M mt C M 3 , Maxwell's equations for the electric and magnetic fields 
E and H and the associated electric displacement D and magnetic flux density B 
are 

(1) curl E(x,t) = —B t (x,t), (Maxwell-Faraday), div B(x,t) = 0, 

(2) cur\H(x,t) = D t (x,t), (Maxwell-Ampere), divD(x,t) = 0. 

Under the assumption of non-dispersive, linear and non-chiral medium, these are 
augmented with the constitutive relations 



Here electric permittivity e and magnetic permeability /i are smooth 3x3 time- 
independent positive matrices. 

The present work consists of two parts. In the first part, we pursue further the 
invariant formulation of Maxwell's equations (QJ)-® for certain anisotropic mate- 
rials, characterized by a scalar wave impedance. This means that the material is 
conformal in the sense that 



where a is a positive scalar function. Physically, the condition of a scalar wave 
impedance is tantamount to a single propagation speed of waves with different po- 
larizations and, without it, the wave propagation gives rise to several generally 
non-Riemannian metrics, see e.g. [23J. Therefore, we consider direct (and inverse) 
problems for the most general subclass of Maxwell's equations which is distinguished 
by the fact that electromagnetic fields with different polarization propagate with 
the same velocity which, of cause, may depend on the propagation direction. The 
scalar wave impedance is encountered in many physical situations. For instance, 
in a curved spacetime with coordinates (x, i) 6 M 3 x K and a "time-independent" 
metric ds 2 = gjk(x)dx' J dx k — dt 2 , Maxwell's equations with scalar permittivity and 
permeability correspond in the coordinate invariant form to Maxwell's equations 
(JIJ-(jH) with scalar wave impedance, see [OJ Sec. 14.1.c] or Sec. 90]. Clearly, 
all isotropic media, i.e. with scalar e and /i have scalar wave impedance. 

The invariant approach leads us to formulate Maxwell's equations on 3-manifolds 
as a first order Dirac type system. From the operator theoretic point of view, this 
formulation is based on an elliptization procedure by extending Maxwell's equations 
to the bundle of exterior differential forms over the manifold. This is a generalization 
of the elliptization of Birman and Solomyak and Picard (see j2j 53j). 



(3) 



D(x, t) = e(x)E(x, t), B(x, t) = fi(x)H(x, t). 



(4) 
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In the second part of the work, we consider the inverse boundary value problem for 
Maxwell's equations with scalar wave impedance. In physical terms, the goal is to 
determine the material parameter tensors e and \x in a bounded domain from field 
observations at the boundary or a part of the boundary of that domain. 

In the invariant approach to Maxwell's equations, the domain M is considered as a 
3-manifold and the vector fields E, H, D and B as differential forms. This alter- 
native formulation has several advantages both from the theoretical and practical 
points of view. First, the invariance of the system and the boundary measurements 
with respect to diffeomorphisms of M that preserve the part of the boundary where 
these measurement are done is essential for the inverse problem. It is possible to 
prove unique identifiability in the invariant formulation and then use this result to 
completely characterize group of transformations between indistinguishable param- 
eters e and \i in the case McK 3 . In particular, when e and \i are scalar functions, 
this result implies uniqueness of the determination of e, /i, and M from data on a 
part of boundary. 

Second, the formulation of electrodynamics in terms of differential forms reflects the 
way in which these fields are actually observed. For instance, flux quantities are 
expressed as 2-forms while field quantities that correspond to forces are naturally 
written as 1-forms. This point of view has been adopted in modern physics, at 
least in the free space, see e.g. [Ej, as well as in applications where the numerical 
treatment of the equations is done using the Whitney elements. An extensive treat- 
ment of this topic can be found in |llU12j . For the original reference concerning the 
Whitney elements, see [Hi] . 

As inverse problems of electrodynamics have a great significance in physics and ap- 
plications, they have been studied starting from the 30's, see e.g. [HHIEE], where 
the one-dimensional case was considered. However, results concerning the multi- 
dimensional inverse problems in electrodynamics are relatively recent. The first 
breakthrough achieved in [SHI UU EEE EI] was based on the use of complex geomet- 
rical optics. These papers were devoted to the identifiability of isotropic material 
parameters e and /i from the fixed-frequency data collected on <9M, namely, the sta- 
tionary admittance map (for the definition of the time-dependent admittance map 
see Definition 11.81 in section 1.3.) Under some mild geometric assumptions it was 
shown there that these data determine isotropic e and /i and also isotropic con- 
ductivity, a, uniquely. These works were based on the ideas previously developed 
in [HI1 IM1 IHj to tackle the scalar Calderon problem, introduced in [13| . Other 
approaches to the isotropic inverse problem for Maxwell's equations work directly 
in the time domain, see [HI [Ell- Regarding the case o = which is considered in 
this paper, the result obtained in proves the identifiability of e and [i from the 
time-dependent data collected on the whole dM in the case when M, considered as 
a Riemannian manifold with metric dl 2 = efj,\dx\ 2 , is simple geodesic. Constructions 



4 



YAROSLAV KURYLEV, MATTI LASSAS AND ERKKI SOMERSALO 



of [8j make it possible to find the product, e/x of unknown parameters e, /x. Moreover, 
the results of (H| are of a local nature making it possible to find this product only 
in some collar neighbourhood of dM. Time-dependent inverse problem for isotropic 
Maxwell's equations was also considered in [TO] which used a time Fourier transform 
to reduce the problem to the one in the frequency domain so that to apply the 
results of |5Dl E]. 

Much less is known in the anisotropic case, where material parameters are matrix 
valued functions. The case of anisotropic e = /x was considered in |H] where it was 
shown that the time-dependent admittance map known on dM makes it possible to 
recover e = /x locally, i.e. in some collar neighbourhood of dM. In spite of very little 
knowledge, it is, however, clear from the study of scalar anisotropic problems that, 
instead of uniqueness, one obtains uniqueness only up to a group of transformations, 
involving proper coordinate changes, see e.g. |H3 EH1 EI3 EH HU US]- A similar 
result for Maxwell's equations was conjectured in [HOI, based on the analysis of the 
linearized inverse problem. Therefore, it is natural to split the study of this problem 
into two steps. First, to formulate and solve the corresponding coordinate-invariant 
inverse problem, i.e., an inverse problem on a manifold. Second, to analyse the 
properties resulting from an embedding of the manifold into R 3 . For a systematic 
development of this approach, see [25J. 

In recent years, inverse problems with data on a part of the boundary have attracted 
much interest, see [H3 EH EB1 EH El] • Part of the motivation come from the physical 
setting when only a part of the boundary is accessible. However, as far as we know, 
there are currently no results on identifiability of the shape of the domain M and/or 
the material parameters e, /x on it from inverse data collected on an arbitrary open 
subset, T C dM. 

A fruitful approach to scalar inverse problems, including those with data on a part 
of the boundary, turned out to be the boundary control method, originated in [5] 
for the isotropic acoustic wave equation. In the anisotropic context, it has been 
developed for the Laplacian on Riemannian manifolds [Jj and for general anisotropic 
self-adjoint [201 EH ) and certain non-selfadjoint inverse problems [33] . 

The current article pursues the study of inverse problems for Maxwell's equations 
significantly further dealing with the global reconstruction of the shape of the domain 
or, more general, 3— manifold M, metric tensor g and scalar wave impedance a, the 
latter two being equivalent to the reconstruction of e and /x. Being based on the 
boundary control, the method developed here combines ideas of the articles [34J 
and [35j with those of [50] and [HI]. What is more, to be able to study anisotropic 
Maxwell's equations, we introduce two essential new ideas. First, we characterize 
the subspaces controlled from the boundary by duality, thus avoiding the difficulties 
arising from the complicated topology of the domains of influence but still providing 
necessary information about the structure of achievable sets, see e.g. Theorem 11.151 
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in section 1.5. This makes our approach much different from that in (HI IH| also 
based on the boundary control method. Indeed, the method of [HI IE] requires local 
controllability in the domains of influence which is no more valid for large times, e.g. 
[0] thus making constructions of [01 IE] inappropriate outside a collar neighbourhood 
of dM. Second, we develop a method of focused waves which enables us to recover 
pointwise values of electromagnetic waves on the manifold and, therefore, reconstruct 
not only the metric g, as in [Hj, but also the impedance a. 

The main results of this paper can be summarized as follows. 

(1) The knowledge of the complete dynamical boundary data over a sufficiently 
large finite period of time determines uniquely the compact manifold en- 
dowed with the travel time metric as well as the scalar wave impedance 
(Theorem I2.1|) . This is valid also when measurements are made on a part 
of the boundary (Theorem 12. 1411 . The necessary time of observation is dou- 
ble of the time required to fill the manifold from the observed part of the 
boundary. 

(2) For the corresponding anisotropic inverse boundary value problem with scalar 
wave impedance for bounded domains in M 3 , the non-uniqueness is com- 
pletely characterized by describing the class of possible transformations be- 
tween material tensors that are indistinguishable from the observed part of 
the boundary (Theorem 12. 18J) . 

(3) For the corresponding isotropic inverse boundary value problem for bounded 
domains in M 3 , the shape of the domain and the material parameters inside it 
are uniquely determined from measurements done on a part of the boundary 
(Theorem EH). 

Some of the results of the paper have been announced in [36| l87|. 

1. Maxwell's equations on a manifold 

This chapter is devoted to Maxwell's equations on a compact oriented 3 manifold 
with boundary. We concentrate on the properties of these equations important for 
the inverse problem considered in Chapter 2. 

We start with the formulation of Maxwell's equations for 1- and 2-forms. These 
equations are augmented to the complete Maxwell system on the full bundle of 
exterior differential forms over a 3-dimensional Riemannian manifold. This allows us 
to define and analyze properties of an elliptic operator related to Maxwell's equations 
and to study the corresponding initial boundary value problem. Crucial results of 
Sections 1.3 and 1.4 are the Blagovestchenskii formula, Theorem ll.91 enabling us to 
evaluate inner products of electromagnetic waves in terms of the admittance map 
Z, also defined in Section 1.3, and the unique continuation result for Maxwell's 
equations with Cauchy data on the lateral boundary. Building on these results, 
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we obtain local and global controllability for electromagnetic waves generated by 
boundary sources and define, in a usual manner, spaces of generalized sources. 

1.1. Invariant definition of Maxwell's equations. To define Maxwell's equa- 
tions invariantly, consider a smooth compact oriented connected Riemannian 3- 
manifold M, dM ^ 0, with a metric g , that we call the background metric. Clearly, 
in physical applications we take Mcl 3 with g being the Euclidean metric. Anal- 
ogously to (II} and ((21), Maxwell's equations on the manifold M are equations of the 
form 

(5) cml E(x,t) = -B t (x,t), div B(x,t) = 0, 

(6) curl H (x, t) = D t (x,t), div D(x, t) = 0. 

Here E, H, D, B e TM, the space of C°°-smooth vector fields on M. They are 
related by the constitutive relations, 

(7) D(x, t) = e(x)E(x, t), B(x, t) = (i(x)H(x, t), 

where e and \i are C°°-smooth positive definite (1, l)-tensor fields on M. We remind 
that, for X G TM, 

(8) (curlX) b = * dX\ divX = -* d* X b . 

Here, d is the exterior differential, b is the fiberwise duality between 1-forms and 
vector fields 

XeTM^X^&M, X\Y) = g (X,Y), 

with VL l M and, generally, Q k M standing for the bundle of differential k— forms on 
M. At last * is the Hodge operator with respect to metric g 0l acting fiberwise, 

* : n k M -> n 3 ~ k M. 

Throughout the paper, we assume that the wave impedance is scalar, i.e., 

(9) fi = a 2 e, 

where a = a(x) is a positive scalar function. This allows us to introduce a new 
metric, g on M. Considering it as a quadratic form on differential 1-forms, we have 

(10) g(X\Y') = -^-g (X\zY'), g * = J-g*^ z = a e = -^ 

As we see later, this metric is responsible for the velocity of electromagnetic wave 
propagation and we call it the travel-time metric. 

Introduce the differential 1— and 2— forms, u 1 , v 1 and u 2 , z/ 2 , 

(11) u l = E\ v 1 = aH\ uj 2 = * B\ v 2 = * aD b . 
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Applying the operator b to Maxwell's equations (5)- (6), we rewrite them in terms of 

UJ 1 , u 1 , uj 2 , v 2 , 

(12) du 1 = -uj 2 , du 2 = 0, 

d-v 1 = -v 2 d-v 2 = 0. 

a a a 

Now, using (fT0|) it is a straightforward matter to check that the constitutive relations 
assume the form, 

(13) v 2 = *uj 1 , to 2 = *u 1 



1 



where * is the Hodge operator with respect to the travel-time metric, g. Eliminating, 
by means of (fl3|) . z/ 1,2 from equations (fT2|. we transform ©-© m t° the system 

(14) w\ = 5 a tu 2 , 5 a uj l = 0, 

(15) u 2 = -dcu 1 , du 2 = 0. 
Here, 5 a : Q k M — > Q 3 ~ k M is the a-codifferential, given by 

(16) 6 a u k = (-if* ad- * u k , 

a 

and * is the Hodge operator in the metric g. These equations are called Maxwell's 
equations for forms on the Riemannian manifold with a scalar wave impedance 
(M,g,a). 

To extend the above equations to the full bundle of exterior differential forms HM = 
Q°M x Q l M x Q 2 M x Q 3 M, we introduce auxiliary forms, uj° e 0°M and u 3 e Q 3 M, 
which vanish in the electromagnetic theory, by 

tO® = SaUJ 1 , U) 3 = —du 2 . 

Since u° = and uj 3 = in electromagnetics, we can modify equations (fT4l) and 
ifnSjl to read 

(17) oj\ = —dujQ + S a u 2 , u 3 = —du 2 , 

(18) u 2 = -dcu 1 + 5 a u 3 , u° = 

or, in the matrix form, 

(19) u t + Muj = 0, 

where to = uj 1 , uj 2 , u 3 ) G OM, and the operator M. (without prescribing its 
domain at this point, i.e., defined as a differential expression) is given as 

/ -5 a \ 

(20) M = I d ° ~ 5a ° 
{ ' d -8 a 

\ d / 
Equations (fl9j) . (|20|) are called the complete Maxwell system. 
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Note that, identifying QM with Q°M ®Q 1 M ®Q 2 M ®tt 3 M , the complete Maxwell 
operator can be written as 

(21) M = d-5 a , 
i.e., it becomes a Dirac type operator on SIM. 
An important property of A4 is that 

M 2 = -diag(A°, Al A 2 , Al) = -A a , 
where the operator A^, acts on fc-forms as 

(22) A* = dS a + 6 a d = A k g + Q k (x, D). 

Here A k is the Laplace-Beltrami operator in the metric g and Q k (x, D) is a first order 
perturbation. Hence, if a; satisfies equation (fT9|) . it satisfies also wave equation 

(23) (d 2 + A a )u = (d t -M)(d t + M)u = 0. 

This formula legitimates the notion of the travel time metric and makes clear that 
in the Maxwell system with scalar impedance, the electromagnetic waves of different 
polarization propagate with the same speed determined by the metric g. 

We end this section with representation of the energy of electric and magnetic fields 
in terms of the corresponding differential forms, setting 

S(E) = ]-[ 9o (eE,E)dV =l [ -co 1 A *a; 1 , 

£{B) = \f 9o ( f iH,H)dV = l [ —u 2 A *u> 2 , 

where dVo is volume form of (M, go) . These formulae serve as a motivation for our 
definition of inner products in the following section. 

1.2. The Maxwell operator. In this section we establish a number of notational 
conventions and definitions concerning the differential forms used in this paper. 

We define the L 2 -inner products for fc-forms in fl k M as 

(to k ,T] k ) L 2 = [ —u k A *r] k , u k , T] k E tt k M, 
Jm " 

and denote by L 2 (Q k M) the completion of Q k M in the corresponding norm. We 
also define 

L 2 (M) = L 2 (tt°M) x L 2 (Q}M) x L 2 (n 2 M) x L 2 (fi 3 M), 
with Sobolev spaces H S (M), Hg(M), sel, given as 

H S (M) = H s (tt°M) x H S (Q}M) x H s (tt 2 M) x H s ({l 3 M), 
H*(M) = H^(n°M) x H^M) x H°(tt 2 M) x H°(n 3 M). 
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Here, H s (Q k M) is the Sobolev space of fc-forms and HQ(Q h M) is the closure in 
H s (VL k M) of the set of fc-forms in f2 fc M, which vanish near dM. 

Clearly, the natural domain of the exterior derivative, d in L 2 (Q k M) is 

H(d,Q k M) = {u k E L 2 (Q k M) | du k E L 2 {Q k+1 M)} , 

and the natural domain of 5 a is 

H(S a ,n k M) = {to k E L 2 (n k M) I 5 a to k E L 2 (fi fc_1 M)} . 

In the sequel, we drop the sub-index a from the codifferential. 

The operators d and 5 are adjoint, i.e. for C^-forms uj k ,r) k+l , 

To extend this formula to less regular forms, let us fix some notations. For uj k E 
Q k M, we define its tangential and normal boundary components on dM as 

tu k = i*u k , nu k = i*(a~ l * u k ), 

respectively, where i* : Q k M — > fl k dM is the pull-back of the natural imbedding 
i : dM — > M. With these notations, Stokes' formula for forms can be written as 

(24) (du k , V k+1 ) L 2 - (u k , 5 V k+1 ) L i = (tu k , n V k+1 ), 

where, for u k E Vt k M and i] k+l E Vt k+1 M, 

(tu k ,nr] k+1 ) = [ tuo k Anr] k+1 . 

JdM 

There are well defined extensions of the boundary trace operators t and n to 
H(d,n k M) and H(5,n k M). The following result is due to Paquet [HI: 

Proposition 1.1. The operators t and n can be extended to continuous surjective 
maps 

t : H{d, Q k M) -> H~ 1/2 (d, Q k dM), n : H (5, Q k+1 M) -> H- 1/2 (d, n 2 ~ k dM), 
where the space H~ 1 / 2 (d,tt k dM) is the space of ' k- forms cu k on dM satisfying 
uj k E H- 1/2 (Q k dM), du k E H- 1/2 (Q k+1 dM). 

Formula (|24j) is instrumental for characterizing the spaces of forms with vanishing 

o o 

boundary data. Introducing H(d, QrM) = Ker (t) and H(S, Vt k+l M) = Ker (n) and 
applying Stokes' formula, one can prove in standard way the following lemma. 

Lemma 1.2. The adjoint of the operator 

d : L 2 {Q k M) D H(d, Q k M) -> L 2 {Q k+1 M) 

is the operator 5 : L 2 (tt k+1 M) D H(5,tt k+1 M) -> L 2 ({l k M) and vice versa. Simi- 
larly, the adjoint of 

5 : L 2 {Q k+1 M) D H{5,Q k+1 M) -> L 2 {Q k M) 
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is the operator d : L 2 (VL k M) D H(d,tt k M) -> L 2 {Q k+1 M). 

When there is no risk of confusion we will write for brevity H(d) = H(d, Q k M) and 
similarly, mutatis mutandis for the other spaces. 

For later reference, we point out that Stokes' formula for the complete Maxwell 
system can be written as 

(25) (77, Mu) L 2 + (Mr], u) L 2 = (tu>, nr?) + (try, no;), 
where to e H with 

(26) H = H(d) x [H(d) n #(5)] x n #(£)] x H(5), 
r] e H 1 (M) , ta> = (to; , to* 1 , tu; 2 ) , nu; = (nu; 3 , nu 2 , nu 1 ), and 

(tu, nr}) = (to; , n^ 1 ) + (to; 1 , mf) + (tcu 2 , nry 3 ). 

With these notations, we give the following definition of the Maxwell operators with 
electric boundary condition. 

Definition 1.3. The Maxwell operator with electric boundary condition, M. e , is an 
operator in L 2 (M), with 

V{M C ) = Ht := H(d) x [H(d] n H(6)) x [H(d) n H(6)) x H(6), 
and M. c uj, uj E V(M. C ) is given by the differential expression (H£J). 

In terms of physics, the electric boundary condition is associated with electrically 
perfectly conducting boundaries, i.e., nxE = 0,n-B = 0, where n is the exterior 
normal vector at the boundary. In terms of differential forms, this means simply 
that tE = tco 1 = and t * B b = tcu 2 = 0. Although not used in the sequel, the 
Maxwell operator with magnetic boundary condition, A4 m , is given by (f271|) with the 
domain 

V(M m ) = H n := H(d) x [H(d) n H(S)} x [H(d) n H(S)} x H(S). 
Consider the intersections of spaces in the definition of V(M. C ) and V(M. m ). Let 

Hl(tt k M) = {u k G H\Q k M) I tu k = 0}, 
HK&M) = {uj k e H\Q k M) I nu k = 0}. 
It is a direct consequence of Gaffney's inequality (see jlZ|) that 

H(d, Q k M) n H(5, Q k M) = Hl{Q k M), 
H{d, Q k M) n H{6, Q k M) = Hl(tt k M). 

The following lemma is a straightforward application of Lemma fL"2l and the classical 
Hodge- Weyl decomposition [57j . 
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Lemma 1.4. The electric Maxwell operator has the following properties: 

(i) The operator M e is skew-adjoint. 

(ii) The operator M e defines an elliptic differential operator in M mt . 

(iii) Ker(A^ e ) = {(0, oj\ u 2 , u 3 ) G H t | dcu 1 = 0, Scu 1 = 0, du 2 = 0, Sco 2 = 
0, 5u 3 = 0}. 

(iv) Ran (M e ) = L 2 {Q°M) x (6H{6, Q 2 M) + dH{d, fi°M)) x 

x (5H(5,tt 3 M) + dH(d,tt 1 M)) x dH(d,tt 2 M). 

By the skew-adjointness, it is possible to define weak solutions to initial boundary- 
value problems needed later. 

1.3. Initial-boundary value problem. In the sequel, we denote the forms w(x, t) 
by cj(t) or to when there is no danger of misunderstanding. 

By the weak solution to the initial boundary value problem 

(27) dtu + Moo = p G Ll oc (R, L 2 W), 

tu} \dMxR = °> ^(0)=u;oGL 2 

we mean the form cj(t) G C(R, L 2 (M)) defined as 

(28) u(t) = U{t)uj + [ U(t- s)p(s)ds, 



where U(t) = exp(— tM e ) is the unitary operator in L 2 generated by M e - Similarly, 
we define weak solutions with initial data at t = T, T G M. Assuming p G C(R, L 2 ), 
the solution has more regularity, uj G C{R, L 2 ) fl C 1 (IR, H'), where H' denotes the 
dual of H. 

We consider also the normal boundary trace nu of weak solutions. Omitting the 
details, we note that nu is defined as a limit of smooth solutions approximating the 
weak one. 

The following result gives a sufficient condition for a weak solution of the complete 
system to be also a solution of Maxwell's equations. 

Lemma 1.5. Assume that the initial dataujQ is of the form uj = (0, Uq, U!q, 0), where 
5(Jq = 0, dujQ = 0, and p = 0. Then the weak solution u(t) of the form $2fy) satisfies 
also Maxwell's equations i.e., uj° = and to 3 = 0. 

Proof: As seen from (|2*3*j) . U)°(t) satisfies the wave equation 

A V + u, t ° f = 0, 

with the Dirichlet boundary condition to; = 0. The initial data for u>° is 
u°{0) = co° = 0, u°(0) = Sco l \ t=Q = 5ul = 0. 
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Hence, u°(t) = for all t G R. 

Similarly, u 3 (t) satisfies the wave equation with the initial data 

w 3 (0) = co° 3 = 0, u 3 (0) = -dco 2 \ t=0 = -dul = 0. 

As for the boundary condition, we observe that 

t5u 3 = tul - tdoo 1 = d t tu 2 - dtuj 1 = 0, 

i.e., the Neumann data for the function *u 3 vanish at dM. Thus, ui 3 (t) = 0. □ 

Assume that u(t) is a smooth solution of the complete system ijTTjl . The complete 
Cauchy data of u(t) consist of 

(tw(x,t),nw(x,t)), (x,t)edMxR. 

The Cauchy data for the solutions u(t) of Maxwell's equations have a particular 
structure. Indeed, by taking the tangential trace of the equation (|T8|) . we obtain 
tu 2 = —dtu 1 . Further, by integrating, 

(29) tu 2 (x,t) =tu 2 {x,0) - [ d(tu\x,t'))dt', x e dM. 

Jo 

Similarly, by taking the normal trace of equation (jTT|) . we find that no^ 1 = dnu 2 , so 
by integrating 

(30) nu) 1 (x,t) = nu 1 (x,0)+ [ d{nu 2 {x, t')) dt', x G dM. 

Jo 

In this work, we consider mainly the case U)(0) = 0, when the lateral Cauchy data 
for the original problem of electrodynamics is simply 

(31) tu = (0,/,- / df(t')dt'), 

Jo 

(32) noo = (0,g, [ dg(t')dt'), 

Jo 

where / and g are functions oft with values in VL l dM. The following theorem implies 
that solutions of Maxwell's equations are solutions of the complete Maxwell system 
and gives sufficient conditions for the converse result. 

Theorem 1.6. If u(t) G C(R, H) nC^M.L 2 ) satisfies the equation 

(33) uj t + Mu = 0, t > 0, 

with u(0) = 0, and tu°(t) = 0, LU 3 (t) = 0, then tu, nu are of the form fHp-fl^). 

Conversely, if tu, nu are of the form ^31\) - K3^I for < t < T, and uo(t) satisfies 
( E^j , with oj(0) = 0, then u(t) is a solution of Maxwell's equations, i.e., oo°(t) = 0, 
u 3 (t) = 0. 
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Proof: The first part of the theorem follows from the above considerations if we 
show that ui{t) is sufficiently regular. For uj 2 G C(R, H(5, fl 2 M)), by Proposition 
mi nu 2 G C(R, H~ l l 2 {Vl l dM)) with drnu 2 G C(R, H'^^dM)). As 5u l t {t) = 
55uj 2 (t) = 0, it holds also that 

mo] G C(R,H- 1/2 {Q 2 dM)), 

implying (JHSJ)- To prove (pTT|) . we use Maxwell duality: Consider the forms 

Then r\ = (r]° , r/ 1 , r/ 2 , r] 3 ) satisfies the complete system dual to the Maxwell system, 
rjt + Air] = 0, where M. is the Maxwell operator with metric g and scalar impedance 
a -1 . Then formula (|32|) for the solution r\ implies (|3T| for to. 

To prove the converse, we observe that the equations, 

(34) d t u°(t) - 5u l {t) = 0, d t Lo l {t) + du°{t) - 5u 2 (t) = 
imply that 

<4 (t) + <Jdw°(t) = 0. 

In addition, uj°{0) = 0, w t °(0) = 0, and from JSO, tcj°(t) = 0. Thus, cj° = for 
< t < T. By the Maxwell duality described earlier, this implies also that uj 3 (t) = 0. 
□ 

The following definition, where R is a right inverse to the mapping t, fixes the 
solutions of the forward problem used in this work. 

Definition 1.7. Let h = (h°,h\h 2 ) G C°°([0, T],CtdM). The solution u{t) of the 
initial boundary value problem 

u t + Muj = 0, t > 0, 

lu(0) = co> G L 2 (M), tw = h, 

is given by 

u(t) = Rh(t)+U(t)u - [ U(t - s){MRh(s) + Rh s {s))ds. 

Jo 

When co>o — and h is a smooth boundary source of the form $31\). 

h = (0,f,- I df(t')dt'), feCfQO^ltfdM), 
Jo 

u(t) is called the solution of Maxwell's equations in M x [0, T] with the boundary 
condition toj 1 = f and the initial condition u(0) = 0. 

To emphasize the dependence of u(t) on / above, we write occasionally 

(35) ^(t)=^(t) = (0,(^) 1 ,(^ / ) 2 ,0). 
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We note that / could be chosen from a wider class, e.g. from if 1/,2 (<9M x [0, T]). 

We use the notation C°°([0, T],Q 1 dM) for the space of C°° functions [0, T] -> Q^dM 
vanishing near t — 0. Theorem 11.61 motivates the following definition. 

Definition 1.8. The admittance map, Z T is defined as 

Z T : C°°([0, T\^dM) -> C^iM^dM), tu l \ aMx[0 ,T] i-> n^ 2 | aM x[o,T], 
where u{t) is the solution of Maxwell's equations in M x [0, T] with uj(0) = 0. 

Note that in the classical terminology of electric and magnetic fields, Z T maps the 
tangential electric field n xE\qmx[o,t\ to the tangential magnetic field nx H\qmx[o,t]- 

The following result, which relates the boundary data and the energy of the electro- 
magnetic field, is crucial for boundary control. It is a version of the Blagovestchenskii 
formula (see jl] for the scalar case). 

Theorem 1.9. (1) For any T > and f, h e C°°([0, 2T], QldM), the knowledge 
of the admittance map Z 2T allows us to evaluate the inner products 

((cufy(t),(uj h y( s )) L 2, 3 = 1,2, o< s ,t<T. 

(2) For anyT>0 and f E C°°([0, T],Vt l (dM)), Z T determines the energy of 
the field uj^ at t = T , defined as 

^(- / ) = ^II(- / ) 1 (t)||! 2 + ^||(^) 2 (t)||! 2 . 

Proof: 1. Let u(t) = oJ f {t) and rj(s) = iv h (s) with F j (s,t) = (uj j (s), rf>(t)) L 2, 
J = 1,2. By (USD, 

(36) (d 2 - d 2 )F^s,t) = {(JMrf®)* ~ {^{s),ri{t)) L , 

= -((d6 + 5d)ui(s),rji(t)) L 2 + (ui(s),(d5 + 6d)rf(t)) L 2 = V(s,t). 

We apply Maxwell's equations (JHJ), and the commutation relations, 

(37) tduo^ = dtu\ nduji = t * *d— * a/ = dt— * = dnu^ , j = 1, 2, 

a a 

where d, in the right side is the exterior derivative on dM. A straightforward 
applications of Stokes' formula (J23J, yields 

b\s,t) = (nu 2 (s),trj\t)) - (tu\s),nr$(t)), 
b 2 (s,t) = (n^^.t^i)) - (tul(s),n V 2 (t)). 

As Z 2T determines b l (s,t) and b 2 (s,t) for t,s <2T and 

(38) F\0,t) = F\s,0) = 0, Fi(0,t)=F/M) = 0, 

the function F^(s,t) can be found from the wave equation (|36|) for s + 1 < 2T. 
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2. Again, by differentiating and using Maxwell's equations and Stokes' formula, we 
obtain 

d t £\ujl) = -{tu\t),mo 2 {t)) = -(/(*), Z T f(t)). 
As S°(ujf) = 0, the energy is readily obtained for t <T. □ 

1.4. Unique continuation results. For further applications to inverse problems, 
in this section we consider the unique continuation of the Holmgren- John type for 
Maxwell's equations. We start with an extension of differential forms outside the 
manifold M. Let T C dM be open and M be an extension of M across T, i.e., 
M C M, T C M int and dM \ V C dM. Let g, a be smooth continuations of g and a 
to M. In this case, we say that the manifold (M, g, a) is an extension of (M, g, a) 
across T. 

Let uj k be a fc-form on M and u k its extension by zero to M. It follows from Stokes' 
formula (2H) that for u k G H(d,Q k M) with ta/| r = 0, we have u k G H{d,Q k M). 

Similarly, if u k G #(<$, fi fc M) and nu k \ r = 0, then u k G #(5, fi fc M). These yield the 
following result. 

Proposition 1.10. Letu(t) G C 1 (M, L 2 )nC(IR, H) with to>|rx[o,T] = and nc<j|rx[o,T] = 
0, be a solution of the complete Maxwell system (CjJ) in M x [0, T]. Le£ u3 6e its 
extension by zero across T C <9M. TTaen u;(t) satisfies the complete Maxwell system 
HjjfyinMx [0,7]. 

We are particularly interested in the solutions of Maxwell's equations. The following 
result extends Proposition II . 101 to this case. 

Lemma 1.11. Let u(t) G C^R, L 2 ) fl C(R, H) be a solution of Maxwell's equations 
OT , <T73j) in M x [0,7], i.e., u°(t) = 0, cj 3 (t) = 0. In addition, let tw 1 | rx [o,ri = 0, 
ncj 2 | rx [o,T] = 0, and w(0) = 0. Then u{t) G C^R, L 2 (M)) n C(R, H(M)) ?s a 
solution of Maxwell's equations in M x [0, T]. 

Proof: The above conditions together with Theorem 11.61 imply that 

to> = (Ojta; 1 , - / dtwW) = 0, ncj = (0, no; 2 , / dnu 2 dt') = 
Jo Jo 

in T x [0, T]. Therefore, by Proposition Qui u(t) satisfies (JE) in Mx [0, T]. Clearly, 
also u°(t) = 0, £ 3 (t) = in M x [0,T], and S3(0) = 0. 

□ 

When we deal with a general solution to Maxwell's equations (|I3 jl - (|I5^ . which may 
not satisfy zero initial conditions, and try to extend them by zero across T, the 
arguments of Lemma fl. Ill fail. Indeed, if u;(0) ^ 0, then (|30|) show that nu; 2 = is 
not sufficient for nu 1 = 0. However, by differentiating with respect to t, the parasite 
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term no; 1 (0) vanishes. This is the motivation why Theorem 11.121 below deals with 
the time derivatives of weak solutions. 

Let, again, T C dM be open and T > 0. Denote by K(T,T) the double cone of 
influence with the base on the slice t = T, 

K(T, T) = {(x,t) e M x [0, 2T] | t(x, T) < T — \T — t\} 

where r(x,y) is the distance function on (M,g) (see Figure 1). 




Figure 1. Left: The double cone of influence. Right: For T large 
enough, the double cone contains a slice {T/2} x M. 

We prove the following unique continuation result for the time derivative of the 
fields. 

Theorem 1.12. Let u(t) be a weak solution of the initial boundary value prob- 
lem ( fg?jJ . Assume that u = (0, 0Jq, Uq, 0), 5ul = 0, oIuJq = 0, and p = 0. If nu 2 = 
on rx]0,2T[, then u t = in K(V, T). 

Proof: When u(t) G C 2 (]0, 2T[, L^nC^O, 2T[, H), then rj(t) = ujj{t) e C x (]0, 2T[, L 2 )n 
C°(]0,2T[, H) also satisfies Maxwell's equations ([T3|) . (|T5|) . Let M be the extension 
of M across T and 77 be the extension of 77 by zero. In follows for (|29| and ijHUj) that 
t^ 2 = — dtcj 1 = and nr? 1 = dnu 2 = in rx]0, 2T[. Therefore, by Proposition lLlO^ 

^(t) G C a (]0, 2T[, L 2 (M)) n C°(]0, 2T[, H(M)), 

is a solution of the complete system and obviously, also a solution of Maxwell's 
equations {TJ), (dHJ in Mx]0,2T[. 

By the unique continuation for sufficiently smooth solutions (see ^H] and Remark 
11.131 below), we have that, for any a > 0, rj = in the double cone 

{(x,t), xeM,te R\t{x,M\M) <T-a- \T-t\}. 
Thus, 77 = in K(T, T). 
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When uj(t) G C°(]0, 2T[, L 2 ) is a weak solution, we use Friedrich's mollifier in t, 

w<j = ip a * u, 4>a(t) = (l/cr)if)(t/a) for a > 0, 

where ip G Cq°([—1, 1]), J tp(s)ds = 1. Then ^(t) satisfies conditions of the theorem 
with rx]0, 2T[ replaced by rx]tr, 2T - ct[. As 

_M^ CT = {-d t yuo a G C°°(]<7, 2T - a[, L 2 (M)), for any j > 0, 

we have in particular that u a (t) G C 2 (](r, 2T - tr[, L 2 ) n C^Qer, 2T - ct[, H). By the 
above, d t oj a (t) = in K a (T,T), where 

K a (T,T) = {(x,t) | t(x,M\M) <T-a- \T-t\}. 
As d t u a (t) — > d t uj(t) in the distribution sense, the result follows. □ 

Remark 1.13. The article [16j, based on results of Tataru [121 EH] deals with scalar e 
and fj,. However, due to the single velocity of the wave propagation, it is, in principle, 
possible to generalize it to the scalar impedance case. Another way to prove the 
desired unique continuation for sufficiently smooth solutions of the equation (|23|) is 
to use the simplified version of Tataru's construction, given in [23 sec. 2.5]. There, 
the unique continuation result is based on local Carleman estimates for the solutions 
of the scalar wave equation, u tt — (lij(x)didjU + Ai(x,D)u = 0, where A 1 (x,D) is 
the first-order differential operator. These estimates utilized a function (f>(x, t) that 
is pseudoconvex with respect to the metric a^, and absorbed the perturbation due 
to Ai(x, D) into the main terms of the Carleman estimates. By (|22|) and (|23|) . 
the operator M 2 is in local coordinates a principally diagonal operator with the 
same second order differential operator, g^d^dj, acting on all components of u)(t). 
As in [2H], one can treat the first-order terms as a perturbation and obtain the 
desired Carleman estimate. In this manner, the constructions in [2H] can be word- 
by-word generalized to the considered case of Maxwell's equations with scalar wave 
impedance. 

Remark 1.14. It is clear from the above arguments that if u(t) is a weak solution 
of the initial boundary value problem $Q and u(t) G C°°(}0, 2T[, L 2 (M)), then 

w(t) G (7 oo (]0,2T[,Z> oo (A'f e )), G C°°(M int x]0,2T[), 

where we used the notation V°°(M e ) = Htvx) ). 

1.5. Controllability results. In this section we derive controllability results for 
Maxwell's equations. We divide these results in local results, i.e., controllability at 
short times and global results, where the time of control is long enough so that the 
controlled electromagnetic waves fill the whole manifold. Both types of results are 
based on the unique continuation of Theorem 11.121 

Let ujf(t), f G C^°(IL|_, Q}dM) be a solution of Maxwell's equations in the sense of 
Definition 11.71 with the initial condition a>^(0) = 0. Let u be the weak solution of 
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(I27j) given by (|2HI) with p = and uj(T) = ujq = (0, Wq, c<j 2 , 0). Similar considerations 
to those in the proof of Theorem II .9[ show that 

(39) (u f (T),u )v = - f (tu f (t),nu(t))dt, 

Jo 

which we will refer to as the control identity. 

1.5.1. Local controllability. In this section, we study differential 1— forms in M gen- 
erated by boundary sources active for short periods of time. Instead of a complete 
characterization of these forms, we show that they form a sufficiently large sub- 
space in L 2 {Q}M). The difficulty that prevents a complete characterization of this 
subspace lies in topology of the domains of influence, which can be very complicated. 

Let T C dM be open and T > 0. The domain of influence, M(T,T), is defined as 

(40) M(r, T) = {x e m\t( X} r) < t}, 

M(r, T)x{T} = K(V, T)n{t = T}. 

Let C°°(]0,T[,fi 1 r) C Cg°(]0, T[, fi^M) consists of forms supported in f x [0, T] 
and 

(41) x(r,T) = c\ L >{{uf)\T) | / e c oo ao,T[,n 1 r)}. 

Furthermore, let 

H(5, M(r, T)) = {u 2 e H(5, n 2 M) | supp (oo 2 ) e M(r, T)}. 

For S C M, we define HQ(Q k S) C //Q(f2 fc M) consisting of fc-forms with support in 
S. 

Theorem 1.15. For any open V C dM and T > 0, 

(42) 5Hq(Q 2 M(T, T)) c X(r, T) c cl L 2 (5H(6, M(T, T))) . 

Proof: The right inclusion being an immediate corollary of (fT4*|) . we concentrate on 
the left one. 

Let e L 2 (fi x M) satisfy 

(43) (u l ,(u f ) l (T)) L2 =0, foraU/6CS o (]0 > T[ j n 1 r). 
By the Hodge decomposition (see |57j) in L 2 (fi 1 M), we have 

(44) ul = ul + 5r,l 

where <£D^ = 0, tu] = and r]^ e H(5, Q 2 M). Thus, (HHJ) is equivalent to 

(45) (^ 2 ,(^) 1 (T)) L2 =0. 

Let uj(t) be a weak solution to (|2T|) with p = and the initial data at t = T given 
by oj(T) = (0, 5r/Q, 0, 0). By the control identity (|39|) . the orthogonality (JUj) and 
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the particular form of the boundary data for solutions of Maxwell's equations (pTTJ) . 
(J32J), we see that 

0= f \tJ{t)^{t))= f (t{u^) 1 {t) 1 mP{t)) = f T (f,nu 2 (t)), 
Jo Jo Jo 

i.e., iiu 2 = on Tx]0,T[. Since 

u(T + t) = (0, u\T - t), -u 2 (T - t), 0), 

also mo 2 = on Tx]T, 2T[. Since 8uj 2 {t) = 0, we see by using Proposition 11.11 that 
i£j 2 G C°(R, H~ l / 2 {Vl l dM)). Hence 

nu 2 = on Tx]0,2T[. 

Therefore, by Theorem ll.121 d t u 2 = in K(T,T), In particular, dSr] 2 , = —d t uj 2 {T) = 
in M(r,T). In other words, if u\ G X(r, T)- 1 , then the term r/ 2 in the decompo- 
sition glj) satisfies cWr/ 2 = in M(T,T). For any i/ 2 G i^(0 2 M(r, T)), we have 
therefore 

(5is 2 ,ujI) l2 = (u 2 ,ddr] 2 ) L2 = 0, 
and thus bv 2 G (X(r,T) ± ) ± . This is equivalent to the left inclusion in (j%2*|) . □ 

Remark 1.16. Later in this work, we deal mainly with the time derivatives of the 
electromagnetic fields. Since u{(t) = uj dt f(t), we see by using 

(46) X(r,T) c cMOj/CT)) 1 I / G (^([O^O^M)} 

and JTJJ), that the inclusions (JUj) remain valid when X(r,T) is replaced with the 
right hand side of (jlfij) . 

1.5.2. Global controllability. This section is devoted to the study of controllability 
results when the control times are large enough so that the waves fill the whole 
manifold. 

For r C DM and T > 0, we define 

(47) y(r, t) = {4(T) | / g Co-Go, n^r)}, 

where VL l Y is the set of 1-forms in Q 1 dM supported on T and abbreviate Y(dM, T) = 
Y{T). Our objective is to characterize Y(T,T) for T large enough. In the following 
theorem, we use the notation 

(48) rad r (M) = maxr(x, T), rad (M) = rad SM (M). 

Theorem 1.17. For open non-empty T C dM and T > T > 2 radr (M), we have 
dh 2 (M)Y(T,T) = Y, where 

(49) Y = {0} x 5H(6) x x {0}. 
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Proof: Let u(t) = ujf(t) be a solution, in the sense of Definition 11.71 of the initial 
boundary-value problem with / G Cq°([0, T ], f^T). Since / = for T > T , we 
have ta; 1 (T) = 0, and consequently, 

u t {T) = -Mu(T) = (0,<5w 2 (T) ^duj^T)^) G {0} x 5H(8) x dH(d) x {0}. 

o 

To prove the converse, we will show that Y(T, T) is dense in {0} x 5H(5) x dH(d) x 

{0}. To this end, let u G {0} x 5H{5) x dH{d) x {0} and u JL Y(T, T), i.e., 

(50) (u; ,u; t (T)) L 2 = {ulu l t {T)) L , + (cu 2 , cu 2 (T)) L2 = 0, 

for any / G C °°([0,T ], M). 

Let a; be the weak solution of the problem 

uJt + M-UJ = 0, tui — 0, u5(T) = u . 

Observe that ojq satisfies Su^ = and du 2 . = 0, so that uj satisfies Maxwell's equa- 
tions. Consider the function F : R -»• R, F(t) = (2(t), u; t (t)) L 2. By Maxwell's 
equations, 

F t (t) = (u(t),u tt (t))v + (u t (t),oj t (t)) h 2 

= —(uj 1 , dduj 1 ) i,2 — (to 2 , d5uj 2 ) l2 + (du , dui 1 ) L 2 + (5cu 2 , 5ui 2 ) L 2, 

and further, by Stokes' formula (J233), 

F t (t) = -(tu 1 ^),^ 1 ^)) - (nu 2 {t),t5uj 2 (t)). 

However, tuj = and 5u 2 = u)\. Thus, 

F t (t) = -(n£ 2 (t),t^(t)> = -(n£ 2 (t),/^)>. 

On the other hand, since u(0) = 0, the orthogonality condition (f50|) implies that 
F(0) = F(T) = 0, i.e., 

(n£ 2 (t), / t (t))di = - / F t (t)di = 0. 
Jo 

Since / G Cg°(]0,T[,fi 1 r) is arbitrary, this implies that nu 2 = in rx]0,T[. Thus, 
by Theorem II . 1 2| £5 tt = in the double cone K(T,T/2). Since T > 2radp(M), this 
double cone contains a cylinder of the form C = Mx]T/2 — s,T/2 + s[ with some 
s > 0. (See Figure 1). 

As uj u satisfies Maxwell's equations with a homogeneous boundary condition tu tt = 
0, this implies that uj u = in M x R. Therefore, cD(t) = c^i + tu> 2 , where ui\ and u; 2 
do not depend on t. Again, by Maxwell's equations, 

L02 — UJ t — M.U)i + tM.U)2- 

Therefore, u 2 = M.uj\ and Aiu> 2 = 0. But then Stokes' formula implies that 

(^2,^2)l 2 = (k>2, Mu)i)x,2 = -(A^W 2 ,^i)l 2 = 0, 
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i.e., u 2 = and Mux = 0. Furthermore, by the choice of u , 
Ul = u{T) =u = (0, -5u 2 , du\ 0) = Mv, 

o 

with v G {0} x H(d) x H{5) x {0}. By a further application of Stokes' formula, 

(ut, ^i) L 2 = (ooi, Mv\2 = -{Mux, z/) L 2 = 0, 
i.e., ljx = and, therefore, Uq — 0. □ 

1.6. Generalized sources. So far, we dealt only with smooth boundary sources 
and the corresponding fields. Later, we need more general sources which are de- 
scribed in this section. 

Let W T be the wave operator, 

W T :C^{]U,T Q [,n x dM)^Y, f h-> uj{{T), 

where T > T > 2rad(M). Let || • ||^- be a quasinorm on the space of boundary 
sources defined via W T , 

(51) 11/11^ =II^ T /IIl 2 . 

By energy conservation, this norm is independent of T > To and by Theorem II .9\ if 
the admittance map Z T is given, we can evaluate ||/||jr for / G Cq°{]0, T [, fl 1 dM). 

Using the standard procedure in PDE-control, e.g. [^31101, there is a Hilbert space 
of generalized boundary sources with the norm defined by (fBTJl . Indeed, we first 
introduce the space JF([0,T ]), 

F([0,T ])=C™(]0,T [,n 1 dM)/ ~, 

where / ~ g iff W T f = W T g, and then complete it with respect to the norm ((ST)) 
to obtain JF([0,T ]). By Theorem I1.17[ W T is an isometry between JF([0,T ]) and 
Y for any T >T > 2rad(M). The elements of JF([0,T ]) are equivalence classes of 
Cauchy sequences (fj)'*L 1 and we denote them by / = (fj)°? =1 . (This is slight abuse 
of notations, as {fj)jL 1 is a representative of the equivalence class /.) To put it in an- 
other way, for any uo G Y, there is a sequence (/j)^ with fj G Cq°{]0, T q [, fl 1 dM), 
defining a generalized source / G JF([0,T ]), and for the corresponding wave 

(52) J t (t) := lim u( j (t), for t > T , 

we have uj{(T) = uq. Since in this work T is considered as a fixed parameter, we 
denote J-"([0,T ]) for brevity as T. 

We say that h G T is a generalized time derivative of / G if for T = T , 

= 0, 

T 



(53) 



lim 
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and write h = D/, or simply h = d t f. We also need spaces with s generalized 
derivatives, T s = V(W), with s G Z + , and J 700 = f] s&+ T s . As in Remark EH if 

/ € T\ 

s 

(54) ufe f| (C s -^[T ,oo[,V(M j e ))nRnn(M e )), 
so that u{{T) G Hf oc (M int ) for T > T . 

We need also the dual of the space V(M s e ). Since C V(M 3 C ), we have (V(M s e ))' C 
H s . Similarly, Hq S C (V(M s c ))'. These facts will be used later to construct 
focusing sequences. 

1.7. Continuation of the boundary data. Theorems 1 1 . 91 and II . 1 71 make possible 
to continue the boundary data, originally given for t < T to larger times t > T, 
when T is large enough, by using essentially the same ideas as in the scalar case, 
[2S1 EH (see also flOj for another continuation method). 

Lemma 1.18. Admittance map Z T , given for T > 2rad(M) ; uniquely determines 
Z l for any t > 0. 

Proof: Let 2e = T-2rad(M). For / G C °°([0,T], Q l dM), Theorem 071 guarantees 
that there is a sequence f n G C^°([e, T], Q}dM) with 

(55) lim uj{ n (T) = ou{(T) in L 2 {Q}M) x L 2 {Q 2 M), 

n— >oo 

or, equivalently, in terms of the energy of the field, 

(56) lim£ T (^") = 0, g n = d t (f-f n ). 

n— >oo 

By using Theorem II .91 one can verify for an arbitrary sequence (/ n )J£Li whether the 
convergence condition (|56|) is valid or not. Moreover, the condition (|56|) is valid 
for some sequence {fn)^=i- Thus, when map Z T is given, one can find a sequence 
(/ n )£Li that satisfies condition (|56|) . 

From the definition (|28|) of a weak solution, (|55|) implies that 

(57) lim n<9 t (u; /n ) 2 |aMx]T,oo[ = n9 i (u; / ) 2 | aMx ] Ti00 [. 

n — >oo 

Let h n (x,t) = f n (x,t + e) G C^°([0,T — e]). Since the function Z T h n determines 
ni9 t (co'^ 1 ) 2 |a Mx ] TiT+e [, we see that Z T determines n(a;^) 2 |9 Mx ] TiT+e [. Iterating this 
procedure, we construct Z l for any t > 0. □ 

In sequel, we need Z T with various values T > 2rad(M). Taking into account 
Lemma fl.!8| we denote simply by Z the admittance map known for all t. 
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Remark 1.19. The controllability results, Theorem 11.151 together with the Bla- 
chovestchenskii formula, Theorem II .91 make it possible to verify from the knowledge 
of Z T whether the condition T > 2rad(M) holds or not. Indeed, T < 2rad(M) if 
and only if for any e > 0, 



This is equivalent to that there are f n G C °°(9Mx]0, T/2[), n = 1, 2, . . . such that 
(u^ n Y(T/2) form a Cauchy sequence in L 2 (M), have norm one and converge to a 
function that is orthogonal to all (w /l ) 1 ((T - e)/2), h E L7 °°(<9Mx]0, (T - e)/2[). 
When Z T is given, this can be verified for all f n and h. 



This chapter is devoted to the inverse problem of electrodynamics. Building on 
properties of Maxwell's equations obtained in Chapter 1, we prove the following 
uniqueness result. 

Theorem 2.1. The boundary DM and the admittance map Z T , T > 2rad(M), 
uniquely determine the manifold M , the travel metric g, and the scalar wave impedance 
a. 

The proof of this theorem consists of several steps. The first is to reconstruct the 
Riemannian manifold (M,g). Having (M,g), we then identify boundary sources 
which generate electromagnetic waves focusing in a fixed point in M at time T > 
2rad(M). These sources are instrumental in reconstructing the impedance a. What 
is more, in section E31 we prove a generalization of main Theorem 12.11 for the case 
where the admittance map is given only a part of the boundary. 

At the end, we return to R 3 to characterize group of transformations of the param- 
eters e and /i leaving the boundary data intact. 

2.1. Reconstruction of the manifold. In this section we determine the manifold 
M and the travel time metric g from the admittance map Z. The idea is to use a 
slicing procedure to control the supports of the waves from the boundary in order 
to determine the set of boundary distance functions. 



M(dM, (T - e)/2) ^ M(dM,T/2). 



2. Inverse problem 



We start by fixing certain notations. Let T < T\ < T 2 satisfy 

T > 2 rad(M) , T x > T + diam(M) , T 2 >2T 1 . 



Let Tj C dM be arbitrary open disjoint sets and < rj < t? < diam(M) be 
arbitrary times, 1 < j < J. We define a set S = S({Tj,r^,T^}j =1 ) C M as an 
intersection of slices, 



(58) 
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(See Figure 2.) Our first goal is to find, in terms of Z, whether the set S contains 
an open ball or not. To this end, we use the following definition. 

Definition 2.2. The set Q = Q({Tj,r^,T^}j =1 ) C consists of generalized 
sources f such that the waves, oo t = lo( , satisfy 



Observe that Maxwell's equations for uj t = w{, f G Q imply that oj tt = (0, b~U)\ , —duj], 0), 
so in particular at t = Xi, we have u)u(T\) = (0, 0, drj 1 , 0), where r/ 1 = —uj]{Ti) has 
the support property supper? 1 ) C S. 



The central tool for reconstruction of the manifold is the following theorem. 
Theorem 2.3. Let S and Q be defined as above. Then: 

(1) If S contains an open ball, then dim(Q) = oo, 

(2) If S does not contain an open ball, then Q = {0}. 

The proof of Theorem 12.31 is given later. 

Theorem 2.4. For any f G J 700 it is possible, given Z, to determine whether f G Q 
or not. 

Proof: Let / = {fk)kLo De a generalized source. By Remark Il.l6| Condition (i) 
of the definition of Q is equivalent to the existence of a sequence h = (/i|)^ , 



(i) ul{T x ) G X{Tj,Tf), for all], 1 < j < J, 

(ii) ^(H) = 0, 

(iii) u tt (T x ) = m M{Tj,t7), for all j , 1 < j < J. 



r 




Figure 2. The set S in the case when J = 2. 




By the linearity of the initial boundary value problem, we have 



\Ku! k ) 1 {T 1 )-{<4*) 1 {T?)\\ = IK^m)!!, 
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where the source is 

9k,e(t) = dt(fk(t) - ht(t + r+ - G Cg°Q0, T x [, tfdM). 

Using Lemma fO^ we can evaluate the norm || (u> 9k >*) l (Ti) || for various (he) and thus 
verify the condition (|59|) . 

In a similar fashion, Condition (ii) is valid for /, if 

lim ||(^) 2 (T 1 ))||=0, 

k— >oo 

which can also be verified via Z by Lemma 1011 

Finally, consider Condition (Hi) for / satisfying Conditions (%) and (ii). Observe 
that 

(d t + M)cJu = in M x R + , 

where u a = oo{ t , and tu u = in <9M x [T , oof. If Condition (Hi) holds, then, by the 
finite propagation speed, 

utt = in Kj = {(x, t) G M x R + \ r(x, Tj) + \t-T 1 \< rf}, 

i.e., u u vanishes in the double cone of influence of Tjx]Ti — r~ ,T X + r~[, for all 
j = 1, . . . , J. Therefore, in each Kj, u t does not depend on time, and, by Condition 
(ii), u 2 = in Kj. Hence, 

(60) nu 2 = Zf = on T, x ] T x - rj , T x + tJ [. 

Conversely, if condition (jfiUl) holds together with Conditions (%) and (ii), then u t 
satisfies 

(d t + M)u t = in M x R + 
with the boundary conditions 

tu] = 0, nuj 2 = in ^ x )T X - tJ , T x + tJ [, 

because Ti — r 3 ~ > T , so that / = in TjX]Ti — t~ , T x + tJ\. Thus, by Theorem 
11.121 ojtt = in Kj and, in particular, Condition (Hi) is valid. As (J60|) is given in 
terms of Z, this completes the proof. □ 

Proof of The,ore.m \&.fk Assume that there is an open ball B C S and let (p G Q 2 M 
with supp(c/>) CC 5. By Theorem 11.171 there is / G T such that 

(61) 4{Ti) = (0,<fy,0,0), 
Clearly, </? G X>°°(jM c ), so that /6 

Let us show that / G Q. Conditions (i)-(ii) are immediate from (j6T|l and Theorem 
11.151 Finally, since 

w£(7i) = -.Mo/ (Ti) = (0, 0, defy, 0), supp(c%) CC B, 
Condition (mj is also valid. This proves the first part of the theorem. 
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To prove the second part, assume that S does not contain an open ball and, however, 
there is / G Q, f ^ 0. Let uj{t) = ujf{t). Then, by Conditions (i)-(ii), 

j 

(62) suppMTO) C f| M(Fj, r+) = S+, 

3=1 

implying, due to u u (Ti) = —Aiu t (Ti), that supp(a>tt(Ti)) C S + . On the other hand, 
by Condition (in), 

u tt {T x ) = in U/ =1 M(Tj, rr) = S~ . 

Thus, supp(a>#(Ti)) C S + \ S~ , which is nowhere dense in M. Since uj tt {Ti) is 
smooth, u u (Ti) = 0, and, therefore, 

(63) duliT,) = -u&T x ) = 0. 
However, by Theorem II .17\ 

^(Ti) = dr] 2 , with ^ 2 G if (5, n 2 M). 

Combining this equation with (|63|) and using tuj\(T\) = 0, we obtain, by Stokes' 
formula pi|l . that 

(^(TO.w^TO)^ = (5r ? 2 ,^ 1 (T 1 )) L2 = (^ 2 , duj](Ti)) L 2 = 0, 

i.e., ct^(Ti) = 0. Also, by Condition (iij, oj 2 {Ti) = 0. These imply that / = 0. 

□ 

We are now ready to construct the set of boundary distance functions, r x , which are 
defined, for any x G M, as continuous functions on dM, 

r x :dM-> K+, r x (z) = t(x, z), z G 9M. 

They define the boundary distance map 1Z : M — > C(9M), 7?.(x) = r x , which is 
continuous and injective, |H21ES1- The set of all boundary distance functions, i.e., 
the image of 

K{M) = {r x G (7(9M) | x G M}, 

can be endowed, in a natural way, with a differentiable structure and a metric tensor 
<7, so that (1Z(M),g) becomes isometric to (M,g), see e.g. jHSEHl. Hence, in order 
to reconstruct M (or more precisely, the isometry type of M), it suffices to determine 
the set 1Z(M). The following result is therefore crucial. 

Theorem 2.5. For any h G C(dM), it is possible, given Z, to determine whether 
h G TZ(M) or not. 

Proof: The proof is based on a discrete approximation process. First, we observe 
that h G TZ(M) if and only if, for any finite subset {zi, . . . , Zj} of dM, there is an 
x G M with 

h(Zj) = T(x,Zj), l<j<J. 
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Denote Tj = h(zj). By the continuity of the distance function, r : M x dM — > R + , 
the above condition is equivalent to the following one: For any e > 0, there are open 
sets Tj C dM, Zj G Tj with diamtTj) < e, such that 



(64) 



wk( f| M(r„ r, + e) \ M(r i; r, - e)) ^ 0. 



On the other hand, by Theorem I2.3| condition (jfiljl is equivalent to 

dim(Q({r j , Tj + e, Tj - e}/ = i)) = oo, 
a condition that can be verified in terms of Z by means of Theorem 12.41 □ 
As a consequence, we obtain the main result of this section. 

Corollary 2.6. The knowledge of the admittance map Z is sufficient for the recon- 
struction of the manifold M and the travel time metric g. 

Having found (M,g), we proceed to the reconstruction of the impedance a. 



oo 

p)p=l 



2.2. Focusing sequences. In this section, we construct sequences of sources, (/ 3 

with the property that (u;/ p ) 2 (T 1 ) = and supp ^(cj/ p ) 1 (T 1 ) \ converging, when 

p — > oo, to a single point in M mt , i.e., the time derivative of the electric field focuses 
to a single point. 

Let y G M mt and 5 y denote the Dirac delta at y in the sense that 
-SJx) A *<f>(x) = <f){y), for any <p G C™(M). 

M « 

Since the Riemannian manifold (M, g) is already found, we can choose Tj P C dM, 
< tJ p < Tjp < diam(M), j — 1, . . . , J(p), so that 

oo 

(65) S p+ i C S p , P| S p — {y}, S p = S{{T j P ,T jp ,T^ p }j=l). 

P =i 

Then, Q p = Q({Tj P , T~ p , t^}^.}) are spaces of boundary sources, which correspond, 
by Definition 12.21 to the sets S p . 

Definition 2.7. For y G M mt . let S p , p = 1,2, ... , be given by A sequence 

(fp)p*Li with f p G Q p , is called a focusing sequence of boundary sources of order s, 
s G Z +; if there is a distribution form A y on M , A y ^ 0, such that 

(66) lim(u;f (T^t?)^ = {A y , V )v, when V G V{M° ). 
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With a slight abuse of notations, we use the same notation for the inner product in 
L 2 and for the distribution duality. We denote a focusing sequences, converging to 
y, by f y = v 

The following theorem characterizes a class of limit distributions that can be pro- 
duced by focusing sequences. This class is large enough for our further goal to solving 
inverse problem. What is more, the sequences from this class can be constructed 
via the admittance map. 

Theorem 2.8. (1) Let y G M mt and (f p )™ 1 be a sequence of boundary sources, 
fp^Qp- Given the admittance map, Z, it is possible to determine, for any s G Z + , 
whether (f p ) is a focusing sequence of order s or not. 

(2) Let f y be a focusing sequence. Then supp(A y ) = {y}. 

(3) For s = 3, the limit distribution A y has the form 

(67) A y =(0,5(X(y)5 y ),0,0), 
where \{y) G A 2 T*M. 

(4) For any y G M mt and X(y) G A 2 T*M, there is a focusing sequence f y , of order 
s = ?,, with (Ay) 1 = S(X(y)5 y ). 

Proof: 1. Take r\ G V(Ail) and decompose it as r] = rji + r/2, where 

rn g V{M s e ) n Y, 7] 2 g v(M s e ) n Y ± . 

As U{t) in (|2*%|) is unitary in V(A4l), by Theorem 11.171 there is a boundary source 

h G T s such that r\\ = oof(Ti). Observe that (uj{ p (Ti),r] 2 )-L2 = 0, so that (f p ) is a 
focusing sequence if and only if there is a limit, 

(68) (A y ,v)v = lim {^{T^^iT^, when h G T s . 

p— >oo 

By Theorem 11.91 the existence of this limit can be verified in terms of Z. 

Conversely, assume that the limit (|68|) does exist for all h G T s . Then, by the 
Principle of Uniform Boundedness, the functionals 



^(7i),77) L 2, peZ 



+ • 



are uniformly bounded in the dual of (T>(Ji4l)) . By the Banach-Alaoglu theorem, 
there is a weak*-convergent subsequence, 

where A y is the sought distribution for which (|66|) is valid. 
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2. Let f y = (f p )^ = i be a focusing sequence. Since f p G Q p , Condition (it) of 
Definition 12.21 implies that A y = (0, Ay, 0,0) and Conditions (i)-(iii), together with 
(ESD, yield 

oo 

(69) supp(dAi) C liminf supp(d(u;f ) 1 (T 1 )) cP|^ = M- 

p=l 

As wf^Ti) G F, 54 = \im p ^ OQ 5{uj{ p )\T 1 ) = 0. Thus, 

(70) supp(A Q Aj) c {y}. 

On the other hand, by Condition (i) of Definition 12.21 

J(p) 

(o^) 1 (T 1 )=0 in M\S p +, S+=f|M(T ip ,r+). 

j'=i 

By the definition (jfiBl) of a focusing sequence, A l y = in M \ S p . As rad(M) < 
diam(M), we can always choose Tj P , t+ so that M \ is non-empty. By the 
unique continuation principle for elliptic systems (see e.g. [20]), it then follows from 
the support property (JTHJ) that supp(A^) C {y}. Since A y is non-zero by assumption, 
supp(Aj) = {y}. 

3. Let s = 3. By part 2. of the theorem, in local coordinates the components 
of A y are finite sums of derivatives of delta-distributions. Since A y G (D(M.l))' c 
H~ 3 (M), it follows that 

3 3 

(71) A\ = ^ (idjdy dx i + Ci5 y dx\ 

i,j=l i=l 

Substituting l|7Tl) into the identity 5^4* = 0, we obtain (|H7j) . 

4. Let ^ p G C^°(5 p ), p = 1, 2, . . ., be 2-forms that converge to A5 y in # x - s (n 2 M). 

By the global control Theorem ll.l7[ there are boundary sources f p such that Loi P {T x ) = 
(0,5ip P ,0,0). Then / y = (f p )^ =1 is a desired focusing sequence. □ 

As y runs over M int , we get a parametrized family of focusing sequences {f y } y( zM int 
which defines the map y 1— > X(y). However, the admittance map does not provide 
a direct access to the values X(y). Although this mapping is unknown, we have the 
following result. 

Lemma 2.9. Given the admittance map Z, it is possible to determine whether the 
map y 1— ► X(y) is a 2-form valued C°° -function in M mt . 

Proof: Let {A y } y£M int be a family of distributions of the form ([67]) corresponding 
to a family {f y } y ^M int of focusing sequences. Assume that y 1— > X(y) is smooth, i.e., 
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A G Vt 2 M mt . Then, for any generalized source h — (hj) G J 700 , we have 

(Ay,^^))^ = {5\5_ y ,^)\T l )) L , = {\5_ y ,d{^) l {T 1 )) L , 

= -(X5 y , (wf 3 ) 2 (T 1 )) L 2. 

By taking the limit j ; — > oo of both sides and using the notation (A, r]) y = *(A A *??) 
for inner product of A, 77 G A k T*M, we arrive at the identity 

(72) (4,,w*(Ti))l> = -<A(y), (wt) 2 (j/,T0) v . 

As (Aj,, cj /i (T 1 )) L 2 = lim^oo lim p ^ 00 (u;/ p ' ! '(T 1 ), : (Ti)) L 2, we can evaluate (f72|) in 
terms of Z by Theorem 11.91 

Conversely, if (A y , u/ 1 ^))^ G C°°(M int ) for any h G J 700 , then A(y) G fi 2 (M int ). 
Indeed, by Theorem II .17| for any <p G f^M, supp(</?) CC M mt , there is a generalized 
boundary source h G with ^(I\) = (0, 0, -cfy>, 0), and, by 

(73) (^,^(Ti)) L 2 = (A(y), d(f(y)) y G C°°(M int ). 

As 93 is arbitrary, condition (f73|) is equivalent to that A(y) is C°°-smooth in M mt . 
□ 

Returning to (f72l) . we conclude that a focusing sequence {/,,} gives rise to a func- 
tional of (u t J ) 2 (T 1 )). It depends only on the value of (u t J ) 2 (T 1 )) at the point y 
and will be called point evaluation functional in the sequel. By the above result 
this functional is determined upto an unknown factor X(y). Hence, by using three 
proper focusing sequences, we can evaluate the 2-form (uj^) 2 at any point in M mt , 
up to a linear transformation. The possibility to control the precise form of this 
transformation is discussed in the next section. 

Lemma 2.10. Let t > T\ and h G T . Given the admittance map Z, it is possible 
to find the 2-forms 

(74) K(y)(u}(y,t)) 2 , y G M int . 

Here K(y) : A 2 T*M -> A 2 T*M is a smooth section of End(A 2 T*M int ) having max- 
imal rank. 

Proof: Let U be a relatively open coordinate patch in M with 2-forms G Vt 2 U , 
k = 1,2,3, linearly independent at any y G U. If {fk(y)}yeu, k — 1,2,3 are three 
families of focusing sequences with the corresponding limiting 2-forms \k(y), we 
define the endomorphism Ku[y) by 

3 

(75) Ku(y)u 2 (y) = ^(A fc (y), u z (y)) y £ k (v), V^U. 

k=i 

As we can evaluate inner products (f72l) by using Theorem II .9| it is possible, for any 
given three families of focusing sequences {fk(y)}yeu, k = 1,2,3 and h, to construct 
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K(y)(u>t(y,t)) 2 for y G U, t > T. Further considerations are based on the result 
that we formulate separately for future reference. 

Proposition 2.11. Let U C M mt be open and G Q 2 Uj k = 1,2,3, linearly 
independent at each y G U '. There are focusing sequences {fk(y)} y eu such that the 
corresponding endomorphism f75| ) is Ku{y) = I y , y G U , the identity in A 2 T*M. 

Proof: Let Xk(y), k — 1, 2, 3 form the dual basis of £k(y), k — 1, 2, 3, 

(^k(y),&(y))y = hi- 
lt is a consequence of Theorem 12.81 that there are focusing sequences giving 
rise to the 2-forms \k(y), which shows the claim. □ 

End of the proof of Lemma V2.1(A By the above proposition, there are, for given lin- 
early independent £k(y), focusing sequences fk(y) so that Kjj(y) is of maximal rank. 
Moreover, since Kuoj 2 {y) can be evaluated for any oj 2 = (u;f) 2 (Ti), the maximality 
of the rank of Kjj(y) can be verified via Z. 

Let Uj, j = 1, . . . , J, be a finite covering of M by coordinate patches and Kj the 
corresponding local endomorphism of the form (j75j) in UjC\M mt . As we can compute 
||Z5|I for all h G J 700 , t > T and x G U, it possible to verify that = IQ(y) for 

y G C/j fl Z7f for all j and I. As by Proposition 12.111 there are families of focusing 
sequences for which this is true, we can construct the desired endomorphism. □ 

2.3. Reconstruction of the wave impedance. So far, we have found the waves 
(u)t) 2 (t), t > T\, up to a linear transformation K which, at this stage, is unknown. 
Since the choice of the focusing sequences is non-unique, we will choose them in such 
a manner that the endomorphism K becomes as simple as possible, i.e., K = CqI, an 
identity up to a constant multiplier. The first step in this direction is to consider the 
polarization of the electric Green's function, defined as the solution of the following 
initial boundary value problem, 

(d t + M)G e (x,y,t) = OinMxt + , 
(76) tG e (x,y,t) = in (x, t) G dM x R + , 

G c (x,y,t)\ t=0 = (0,5(A^),0,0). 

Sometimes, we denote G e (x, y, t) = G e (x, y, t; A) to indicate the source A G A 2 T*yM. 
Assume that h = h y is a focusing sequence that produces a wave focusing at y, the 
corresponding 2-form being A. Since the boundary sources are off when t > T\, we 
must have 



(77) 



G e (x,y,t) = u*(x,t + T 1 ). 
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On the other hand, by Lemma l2.10[ we can calculate the 2-forms K(x)(uj^) 2 (x, t+Ti) 
for x G M mt and t > 0. Hence, we know the electric Green's function up to a linear 
transformation. 

Let us denote by $ = $ A (x, y,t) the standard Green's 2-form, satisfying 
{d 2 t + A 2 a )$(x,y,t) = 0inMxl + , 

(78) $(x,y,t)\aMxVL+ = 0, 

(79) y, 0) = 0, $ t (x, y, 0) = \(y)5 y (x), 

where A(y) G A 2 T*M and the boundary condition in l|7nj) means that all three 
components of $ vanish on dM x R + . 

Let G e = G e (x,y,t) be defined as 

(80) G c = (d t - M)(0, 5$, 0, 0) = (0, d t 5$, -d5Q, 0). 

As (d 2 + A a ) = (9 t + *M)(<9i — .M), G e satisfies the complete Maxwell system and, 
by (f79|) . the initial condition in (fTBj). By the unit propagation speed, G e = near 
9Mx]0,r(?/,9M)[, satisfying the boundary condition in (JTBJ) . Thus, G e (x,y,t) = 
G e (x,y,t)fmt<T(y,dM). 

To further the study of G e , we formulate the following result proved in the Appendix. 

Lemma 2.12. For every y G M mt t/iere is an open neighborhood U C M int of y, 
a positive t y and a mapping Q y (x) that is smooth with respect to x G U , where 
Q y (x) : A 2 T*M -> A 2 T Z *M «s bijective, such that 

(81) $(x,y,t) = Q y (x)\S(t 2 -T 2 (x,y)) + r(x,y,t), {x, t) G C/x]0, t y [. 

Moreover, with some smooth Q? y {x) : A 2 T*M -> A 2 T*M, p = 1,2 and C^-smooth 
2-form r(x,y,t), the remainder can be written as 

(82) r(x, y, t) = £ W A(* 2 - r 2 (a;, y))^ + f(x, y, t). 

p=l,2 

By (|80| . it follows from (|8T|) that, for sufficiently small t, 

(83) G c = (0,Gl,G 2 e ,0) + ri , 
where 

G£ = -2t * (dr 2 A *Q,,A)5 (2) (t 2 - r 2 ), G 2 = dr 2 A *(dr 2 A *Q,,A)5 (2) (£ 2 - r 2 ), 

and T\ is linear combination of a bounded function, the delta distribution on dB y {t) 
and its first derivative, B y (t) being the ball of radius t centered in y. Using Lemma 
12.101 and (|77jl together with (|83|) . we obtain the following result. 
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Lemma 2.13. Given the admittance map Z, it is possible to find the distribution 
2- form 

K(x)G 2 c (x, y, t) = K(x)(u}) 2 (x, t + 71), 

where K G End(f2 2 M int ) and t > 0. Moreover, the leading singularity of this form 
when < t < t y determine the 2-form 

(84) K(x) (dr 2 (x, y) A *(dr 2 (x, y) A *Q y (x)X)) , x G dB y {t). 

The linear transformation K(x) of Lemma EHH depends on fk(x), £k(x), k = 1,2, 3. 
Our next goal is to formulate conditions, in terms of Z, on to make K isotropic, 
i.e., 

(85) K{x) = c(x)I, c G C°°(M int ), c(x) + 0. 
To this end, observe that for A G A 2 T*M, 

(86) t By {t) {dr 2 A *(c/r 2 A *Q y X)) = 0, 

where ts (<)C<;* is the tangential component of u k on dB y {t). Physically, condition 
(jKKl) corresponds to the orthogonality of the polarization of the magnetic flux density 
and the direction of wave propagation. (See Figure 3.) If K is isotropic, we have 

(87) t By{t) [K{dr 2 A *(dr 2 A *Q y X))j = 0. 




Figure 3. Vector v is the right polarization of the electromagnetic 
wave in the plane M x {t}. The reconstructed polarization w has 
wrong direction, if the transformation matrix K{x) is not isotropic. 

Conversely, we show that condition ((87]) for all y G M mt and t <t y guarantees that K 
is isotropic. What is more, the condition (|87]) is verifiable from the knowledge of Z. 
Indeed, for \{y) G A 2 T*M and t = r(x, y), JHZ|) means that K(dr 2 A*(dr 2 A *Q y X)) 
is normal to T x dB y {t) C T X M, i.e, for vectors X, Y G T x dB y (t), we have 

K(x)(dr 2 (x,y) A *(dr 2 (x,y) A *Q y (x)X(y)))(X,Y) = 0. 

Observe that when X(y) runs through A 2 T*M 7 then *(Q 2/ (x)A(y)) runs through 
T*M. Now we may vary y and t with fixed x such that t(x, y) = t, making T x dB y (t) 
run through the Grassmannian manifold G 3;2 {T X M). Transformation K(x) is kept 
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invariant in this variation. Hence, we deduce that K(x) keeps any 2-dimensional 
subspace of A 2 T*M invariant, so it must be isotropic as claimed. 

Assume that the the focusing sequences used for the point evaluation functionals 
are chosen so that K(x) = c(x)I. For any generalized source / G J 7 , we may thus 
evaluate 

(u!) 2 (x,T 1 ) = c(x)(lu!) 2 (x,T 1 ), 
with yet unknown c(x). Since w* satisfies Maxwell's equations, we have 

d(u{) 2 = dc A {J t f + cd(J t f = dcA (oof) 2 . 
The global control Theorem II .171 thus asserts that c(x) = c is equivalent to 

(88) d{u(f{x, Ti) = 0, for all / G F, 

a condition that is verifiable from the knowledge of Z. Hence, the focusing sequences 
used for point evaluation can be chosen such that c(x) = Co ^ 0. 

To proceed with reconstruction of a, consider the inner product, 

(89) / (£/) 2 (*,T\) A *(2?) 2 (a ?J T 1 ) = c% [ (u;/) 2 ^, 1\) A *{u}) 2 {x,T x ), 

JM JM 

which can be found via Z. On the other hand, by Theorem II .9[ Z determines the 
energy inner product, 

2 J M a(x) 

By choosing the boundary source h = hj such that h = (hj) ( *L 1 is a focusing sequence 
and by comparing the above inner products at the limit j — > oo, we recover the value 
Cq«(x) at any point x G M. 

Finally, we notice e.g. by considering the energy integrals that the admittance map 
has the scaling property Z^ Mg c 2 a ^ = c~ 2 Z^M,g,a)y with evident notations. Therefore, 
given Z and (g, c^a) already reconstructed, it is also possible to determine Co and 
hence a. This completes the proof of Theorem l2.il □ 

2.4. Data given on the part of the boundary. In this section, we generalize 
the proof of Theorem 12. II for the case when the data is given on a non-empty open 
subset T C dM. In this case, instead of the complete admittance map Z T we are 
given the local admittance map Z^ , defined by 

zif = z T f\ rxnn , f g c^acn^r), 

where is the space of 1-forms / G fl 1 dM supported on V. Denote Z = Z T with 
T = oo and recall that radr(M) is the geodesic radius of M with respect to T, see 

(Hi. 
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Theorem 2.14. Given T, the local admittance map Z-f , T > 2radr(M) ; uniquely 
determines the manifold M, the metric g, and the scalar wave impedance a. 

Proof. Here we use notations of section II .61 By Theorem II .171 we have that the set 
T v = C^°(\0,T o [, Q^Jy with T > T > 2 rad r (M) is a dense subset of T. Thus we 
can identify with T . This makes it possible to use, when the data is given on a 
part of the boundary, all the results about generalized sources obtained in section 
II .61 for the whole boundary. In particular, we can define sets C J~r that can be 
identified with T s . 

Exactly as in section 11.41 we can show that the local admittance map Z^ , T > 
2 radr(M), determine the map Z\ for all t > 0, i.e., the map Z?. 

Our first aim is to reconstruct (M,g) near T. For this, let 9 T : V — > R be 
9 r (z) = sup{s > : T(y ZtU (s),T) = s} 

and 

M r = {y z ,„(s) G M : zeT, 0<s< 6 r {z)}, 

where j Z)V (s) is the geodesic starting from z G dM in the normal direction. 

Lemma 2.15. Given T, the local admittance map Z-p determines the function 6 r : 
r — > K., the Riemannian manifold (M-p,g), and the wave impedance a on Mp. 

Proof: Let z G T. Using notations of section l2~T1 we see that s < r (z) if and only if 
S = M(T Z , s) \ M(r, s — e) has non-empty interior for all open Y z C V containing 
z and e > 0. In turn, S = (M(T Z , s) \ M(T Z , 0)) n (M(T, s) \ M(T, s - e)) is of the 
form (|58|) . Thus, using Theorem 12.41 with Z r instead of Z, we can find out whether 
S has non-empty interior or not. Thus we can find r (z). 

At this stage, we can proceed in a similar manner to the proof of Theorem 12.51 to 
find functions Rr(M r ) = {r x \ r : x G M r }. To do this, we just use the procedure 
presented in the proof of Theorem l2.5l but consider only functions h G C(T) that have 
a unique global minimum, say z G T, with h(z ) < r (z ). After construction of 
this set, we see as in [23 Sect. 4.4] that the set Rp(Mr) determines the Riemannian 
manifold (M r , g). 

Reconstruction of a in M r follows the same route as with data given on the whole 
boundary by restricting our attention to focusing sequences corresponding to points 
y G M r and using identification of jFp with T s . □ 

In the next step we will show that we can find the admittance map on the boundary 
of an arbitrary ball B C M r . We will denote by Zq B the local admittance map 
defined by using the manifold M\B instead of M and dB instead of Y. For similar 
arguments in the scalar case, see [26] . 
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Proposition 2.16. Given (Mr,g,a) and the map Zr for (M,g,a), we can find the 
local admittance map Zqb for (M \ B, g, a). 

Proof. First we observe that (Mr, g, a) and Zr determine the values of the electric 
Green's function G 2 (x, y, t; A) for any x, y G M r , t > 0, and A G A 2 T*M. This result 
is proven in section l2~3l in the case when the admittance map is given on the whole 
boundary and the proof can be directly extended to the considered case. 

Consider now the initial boundary value problem 

(90) d tV + M V = K, lxiMxR + , trj\ dMxR+ = 0, tj(0) = 0, 

where k = (0,8(3,0,0) with (3 G C°°(IR+, £l 2 B) and denote its solution by rj = rjp — 
(0, rj 1 , rj 2 , 0). Writing rf in terms of the electric Green's function we obtain 

(91) r, 2 {x,t) = [ [ G 2 e (x,y,t-t';^^-)dV g (y)dt', 

Jr+ Jb <x\y) 

where dV g is the Riemannian volume measure on (M,g). Using equation (|9"UJl we 
also find 

(92) n\x,t) = [ {5r] 2 {x,t') + 6/3{x,t'))dt'. 

Jo 

We continue the proof with the following lemma; 

Lemma 2.17. Let uj = (0, uj 1 , uj 2 , 0) G (^(R+.fiM) be a solution of Maxwell's 
equations j tiTj] and j fj^jj in (M \ B) x K + which satisfies the electric boundary 
condition too = on dM x M + and initial condition uj(0) = 0. Then there is 

(3 G C°°(IR+,^ 2 S) such that the solution rip of initial boundary value problem fiJTA) 
coincides with uj tt in (M \ B) x R + . 

Proof. Let S = (0, S 1 , cD 2 , 0) G (^(IR+^M) be an arbitrary smooth continuation 
of uj into Bxl + . Let 

p = (0, p 1 , p 2 , 0), p 1 = -8dw l , p 2 = -dbZ) 2 . 

Then p = uj tt in (M \ B) x R +) and p satisfies Maxwell's equations 

p* — 5p 2 = (Ma 1 , a 1 = — uj t + 8ui 2 , 
p 2 + dp 1 = d5a 2 , a 2 = —uj 2 — duo 1 , 

in M x R + . Then rj = (0, p 1 — 8a 2 , p 2 , 0) satisfies the initial boundary value problem 
dSU with /3 = da 1 —a 2 supported inBxl + . In particular, uj u = r)p in (M\B) xl + . 
□ 

To complete the proof of Proposition ^. 16l we start with an arbitrary (3 G C°°(K + , Vt 2 B) 
and find, using formulae (|9T| . (|92|) . the wave ^(x, t) for x G M r . Let u;(x, t) be now 
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defined as 

(93) u°(x,t) = 0, u l (x,t)= f [ r]l(x,t")dt"dt', 

Jo Jo 

t r t' 

"\ Jj-ii Jj-i 



JO 



u 2 (x,t) = / / rfp(x,t")dt"dt', uj 6 (x,t)=0. 



Then u(t) is a solution of the initial-boundary value problem, 

u t + Mu = in{M\B)xR + , u{0) = 0, 

tk>|aMxR+ = 0, tuj\Q BxK+ = (0, fp, — / dfp), 

Jo 



where 

(94) fp= f [ tr]p(x,t")dt"dt' ecr^L+^dB). 

Jo Jo 

Using again formulae (|9TJ) . (|92|) . we see that (M r ,g, a), Z r determine the map 

(95) fp ^ nu 2 \ dBxR+ = [ [ nr)l(x,t")dt"dt' Grfl+.Q^) 

Jo Jo 

for any fp of form (jHH). As according to Lemma l2.171 the map (3 — > fp is a surjective 

map from C°°(^+, & 2 B) onto C 00 ^, Vt l dB), the map (gH) determines Z T . This 
proves Proposition 12.161 □ 

Having found Zqb we construct the Riemannian manifold Mqb C M \ B, metric 
g on Mq B and impedance a on Mg B . Here Mq B is defined in a similar way as M r 
changing M to M \ B and T to 95. Combining this with the previous results, we 
find the part M r U Mq B of M as well as the metric g and the wave impedance a 
on it. Iterating this procedure, we reconstruct, in finite number of steps, the whole 
manifold (M,g, a). For detail, see [23 Sect 4.4.9]. This proves Theorem f2. 141 □ 



2.5. Inverse problem for Maxwell equations in IR 3 . In this section, the unique- 
ness results for Maxwell's equations on a manifold are used to characterize the non- 
uniqueness of inverse problems for Maxwell's equations (EQl-Q in a bounded domain 
of M C K 3 with the Euclidean metric (<?o)r) = 

Let Mj C M 3 , j = 1, 2, be two bounded smooth closed domains with a common part 
T of their boundaries, T C dM 1 fl dM 2 . Let ej and /Xj, j = 1, 2 be permittivity and 
permeability matrices in Mj, respectively, with [ij = a 2 ej, with <x,- > being the 
corresponding scalar impedances. Assume that the local admittance maps Z r j for 
(Mj, €j, Hj) coincide. By Theorem l2.14[ both (Mi, e\, fix) and (M 2 , ei, /i 2 ) correspond 
to the same abstract manifold (M, g, a) which is uniquely determined by Zrj with 
the part T corresponding to a part T C <9M. This implies that there are embeddings 
Fj ■ : M — > Mj C M 3 of the manifold M in the Euclidean space such that the 
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metric tensors and the wave impedances satisfy g = (Fj)*gj and 5 = (Fj)*acj and 
Fi|p = F 2 \f. Recall that gj are determined by expression (|TTH) with Sj and jij in 
place of e, \i. Embeddings Fj induce a diffeomorphism 

(96) $ = F 2 o Ff 1 : Mi -> M 2 , $| r = id. 

Consider two vector fields X x and Y"i in Mi, and denote X 2 = D&Xi, Y 2 = D<&Y\. 
The electric energy inner product for the corresponding 1-forms u\ rj 1 G VL l M is 
invariant, i.e., we have 

/ g (X 1 ,e 1 Y 1 )dV = f 4a; 1 A V = / So(X 2 , e 2 F 2 )^ . 

On the other hand, as X 2 = X 1 and Y 2 = D^Y X , 

[ g (X 2 ,e 2 Y 2 )dV = [ g (X 1 ,$*e 2 Y 1 )dV , 

where 

(97) <Te 2 = 5 -^L_( jD $) T (e 2 o $) 

Since X\ and Y"i are arbitrary, we must have e% = $*e 2 . Similar reasoning shows 
that iii = $*/i 2 - 

Thus we have proven the following result. 

Theorem 2.18. Let Mi, M 2 C R 3 be bounded smooth domains and T C 9Mi fl<9M 2 
oe open and non-empty. Let 2r,i <™d ^r,2 fre t/ie local admittance maps correspond- 
ing to (Mi,ei,/ii) anc? (M 2 ,e 2 ,/i 2 ), respectively. Then .2^,1 = -2^,2 on/?/ «/ 
i/iere is a diffeomorphism $ : Mi — > M 2 , $L = id and €\ = $*e 2 , [i\ = $*yU 2 . 

Remark 2.19. It follows from (J9"TJ) that e and /i do not transform like tensors 
of type (1,1). This is due to the special role played by the underlying Euclidean 
metric g$ = 8^ , which is not changed by diffeomorphisms $. These transformations 
were observed also in the study of the Calderon inverse conductivity problem. It is 
shown in [HHl that, for Del 2 , boundary measurements determine the anisotropic 
conductivity up to same group of transformations as described in Theorem 12.181 
For n > 3, a similar result is conjectured, based on the analysis of the linearized 
inverse problem, see jHDl- The Calderon problem is closely related to the inverse 
problem for Maxwell's equations, as the low-frequency limit of Z is related to the 
Dirichlet-to-Neumann map for the conductivity equation [42J. 

When e and \i are isotropic, we obtain the following uniqueness result. 

Theorem 2.20. Let M C M 3 be a bounded smooth domain, T C dM be open and 
non-empty, e and /i be smooth positive functions on M and Z? be a local admittance 
map for (M, £, ji). Then V and Z-p determine (M, e, ji) uniquely. 
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Note that the knowledge of M is not a priori assumed in the above theorem. 

Proof. Assume that for (Mi,ei, pi) and (M2, e 2 ,// 2 ) such that T C dM\ fl <9M 2 we 
have Zr,i = 2r,2- Then there is a diffeomorphism $ : Mi — » M 2 satisfying <3>| = id 
and t\ = $*e2, pi = $*/i2- Since ei and 62 are isotropic, it follows from the Liouville 
theorem that $ is conformal. Since $L = id, it follows that $ is identity. □ 

2.6. Outlook. There are several direction to which the present work can be ex- 
tended. 

1. A natural inverse problem is the inverse boundary spectral problem for the 
electric Maxwell operator A4 e . The problem is to determine the metric g and wave 
impedance a, or, in the other words, e and p from the non-zero eigenvalues Xj of M. e 
and the normal components of the corresponding eigenforms on dM. This problem 
was studied in, e.g. |H], for the scalar Maxwell's equations. For the considered 
anisotropic case, this requires significant modifications of the method developed in 
this paper and will be published elsewhere. 

2. With the uniqueness of the inverse problem in hand, the next issue is to study 
stability of the inverse problem and develop stable reconstruction algorithms. A 
general approach to these questions, in the scalar case, is introduced in (2^j, in terms 
of certain geometrical a priori bounds on (M,g), with sharp results on conditional 
stability in pp. Adding a priori analytical bounds on a, we intend to analyse these 
questions for anisotropic Maxwell's equations. 

Appendix: The WKB approximation. Denote by $(x,y,£) = $\(x,y,t) the 
Green's 2-form, i.e., the solution of 

(98) (d 2 + A 2 a )$ x (x) = in M x R + , 

®\{x)\t=a = 0, d t <£> x (x)\t=o = X5 y (x), <S>\(x)\ dMxR+ = 0, 

where A G A 2 T*M. Let B y (p), p < r(y,dM) be a domain of normal coordinates 
centered at y, so that 

(99) g»(0) = 5^, d k gv{0) = 0. 

Rewriting equations (J98)) . componentwise, in these coordinates and using the unit 
propagation speed, we can, instead of (|98j). consider the fundamental solution, 

$(x,y,t), t < p, 

(100) {{d 2 t - g ij didj)I + S^ + C}$ = 0, inMx]0,p[, 
$|t=o = 0, $$|t=o = 15{x), 

where I is the 3x3 identity matrix and B l (x), C(x) are smooth 3x3 matrices. 
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Following [El IS], which deal with the scalar case, we search for the solution to 
in the WKB form: 

(101) t) « G (x) 5(t 2 -r 2 ) + J2 Gt(x) {t 2 - r 2 y- l /{i - 1)!, 

where r(x,y) = \x\. Substitution of (|1()1|) into equation (IK)Oj) gives rise to the 
recurrent system of transport equations. The principal one is the equation for Go, 

dC 

At^(tx) + {(g ij (rx) d^r 2 - 6) / + BHrx) ^r 2 } G (tx) = 0, 

where x — x/r. To satisfy initial conditions in ()100|) . we require that G (0) = 
(2n)- 1 I. By ® , g^d^T 2 }^ -6 = 0. Also, d { T 2 \ x=0 = 0. Therefore, 

-!- {{g ij {rx) d t d 3 r 2 -6)1 + B\rx) d.r 2 } 
4r J 

is a smooth function of (r,x), so that Gq(x) is a smooth 3x3 matrix- function of 
(r, x), for r > 0. Moreover, it can be shown that Go(x) is also smooth at x = 0. 

For G^, £ > 1, we obtain transport equations 

4r^ + {{U - 6 + g ij (x)d i d j T 2 (x)) I + ^(s^r 2 } 

with G^(0) = 0. If we write G? = GoFi, we obtain for Ff the equations 

(102) 4r^ + 4£F, = G" 1 [^'9^1 - B% - C] G^, F t (0) = 0. 
Solving the above equation, we find 

(103) F i (x) = ]r- e f G \sx) {[g^d^I - B%- C] G^} (sx) s^ds, 

4 Jo 

which are smooth functions of x. As (llOOj) is a hyperbolic system, the right side of 
(llOlj) is the asymptotics, with respect to smoothness, of $(x, y,t), when t < p. 

Clearly, the asymptotic expansion (IK)ljl implies decomposition (|%T1) . (j%2~| . 
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